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Abstract 

This is the fifth paper in the series outlining an algorithm to consis- 
tently quantize four-dimensional gravity. We derive the pure Kodama 
state by path integration, in analogy to the no-boundary proposal for 
constructing quantum gravitational wavefunctions, checking at each 
stage of the process the equivalence of the canonical and path integral 
approaches. A family of additional pure Kodama states is identified 
via the canonical approach and a criterion for their suitability as a ba- 
sis of states is examined. We provide an interpretation for the problem 
of time within the context of generalized Kodama states and propose a 
possible method of resolution. We also develop different techniques for 
solving the Gauss' law constraints at the kinematical level, in prepa- 
ration for future work in this series. 



1 



1 Introduction: The no-boundary approach to quan- 
tum wavefunctions 



The no-boundary proposal is a path integral representation for the wave- 
function of the universe, based upon the transition amplitude of the 3-metric 
hij = hij {x, t) between an initial configuration on a 3-surface Eq and a final 
3-surface Ey, including a matter field (f) = (f){x,t), given by [8]. 

^HH[hij{^T)A{^T)]= j dg^,D^ei,v[-SEH{a^..A)/h] (1) 

where the gravitational portion of the path integral is over all 4-metrics 
g^ii, throughout the interior of a compact 4-manifold M with the given 3- 
metric and matter fields {hij{'ET),4>{^T)) and (/lij(Eo), ^(Eq)) induced on 
the spatial 3-boundariy Et and Eq. The action for the gravitational system 
is the Einstein-Hilbert action, given by 



Ieh = [ (fx^i^^^R - 2A) + S[(t>i f d^xVhivK, (2) 

iDvrG Jm Jt. 

where g^^ is the 4-metric of spacetime and Kij is the extrinsic curvature of 
a 3-surface E of intrinsic curvature ^^^R. To find the Hamiltonian for the 
system (and thus the energy eigenstates) , one performs a 3-1-1 decomposition 
of the Einsten-Hilbert action (refs[8],[5]) 

S[g] = [ {T^'^kj +Tr<P -NH- N'Hi) (3) 

JM 

where H and Hi are the classical Hamiltonian and diffcomorphism con- 
straints corresponding to and N'^^ the lapse function and shift vector, 
respectively, given by 

H = -GijknT''T^^^ -Vh^'^^R + ikVh ■ Hi = 'iTij\j, (4) 
where Gijki, the metric on superspace, is given by 

Gijki = l/2Vh{hikhji + hiihjk - hijhu) (5) 

and the momentum conjugate to the induced 3-metric hij on E, namely tt*-', 
is given by vr*-^ = G^^'^^Kki- 

Since the gravitational Hamiltonian is a linear combination of first-class 
constraints, it must be the zero of energy for the system, which corresponds 
to an exact cancellation of the gravitational against the matter energy, at 
least for compact manifolds. This is enforced by the requirement that the 
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wavefunction be invariant with respect to variations of the Lagrange multi- 
phers. Define ^'0(5, ^) = ^pathig^^P]- Then we must have 



Sj^^^pathig, 4>] = sj;^%ath[g, 4>] = o Vx. (6) 

Both the path integral and canonical approaches to quantization con- 
tain some elements of ambiguity, namely operator ordering ambiguities in 
the former and ambiguities in the path integral measure in the latter, but 
nevertheless it has been established that the two approaches must, at least 
formally, be equivalent to one another ([8], [9], [10]). 

For some simplified cases, reasonable equivalence has been shown to exist 
between the two approaches in minisuperspace models for certain restrictions 
upon boundary conditions, operator ordering and convergence ([8], [9], [5]). 
Still, one would like also to be able to address these issues within the context 
of the full theory, unobscured by any simplifications due to minisuperspace. 

The path integral measure can be decomposed into the following form 



DgD(P= dN{x)dN\x)dhki{x)d(p{x) = DNDNiDhDcp. (7) 

The purpose of this decomposition is to separate the gauge (A'', iV*) from 
the physical {hij , cj)) degrees of freedom in view of the fact that they both 
originated from the same total phase space {gnv^4>), and to make it more 
physically clear the sequence of path integration along these variables. 

The invariancc of the wavefunctional with respect to variations in the 
Lagrange multipliers corresponds to the implementation of the constraints 

^cnn[g, 4>] Hi{x)^can = Vx (8) 



5N\ 

which is the diffeomorphism constraint, independently valid at point a; of M 
and 

-*can[5, (P] H{x)^can = Vx (9) 



SN{x) 

which is the Hamiltonian constraint, also independently valid at each point 
x of M. 

We will define an analog to the no-boundary proposal for determining 
quantum gravitational wavefunctions in Ashtekar variables. The argument 
of the wavefunction must include, as a minimum, the values of the fields 
(^,«(So),</.(So)) and 

{Af{T,T), </)(Sr)) fixed on the respective 3-boundarics. One obervation, as in 
the no-boundary proposal [8], is that the additional degrees of freedom path- 
integrated within the interior of M beyond those defined on the boundaries 
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Et and Sq are pure gauge degrees of freedom. In the case of the Ashtekar 
variables one gauge degree of freedom is the time component of the connec- 
tion Aq = 0'^. Note that the densitized lapse N_ and the shift A^*, strictly 
speaking, are not part of the original phase space of the Ashtekar variables 
in the same manner that they are for the metric variables, in the sense that 
they do not combine to form a covariant object such as the 4-metric gfj_y. In 
order to impose all corresponding constraints within the interior of M, the 
phase space must be 'enlarged' to include these 'gauge' degrees of freedom. 
So, in Ashtekar variables we might expect a relation of the form (for pure 
gravity without matter, for simplicity), schematically. 



to yield the wavefunction. A relevant question becomes what version of 
the constraints is imposed by the delta functionals arising from the path 
integral: the classical or the quantum version. 

Since the path integral involves c-numbers and not operators one may 
attempt to attribute the 'knowledge' of ordering ambiguities arising from 
the canonical, operator approach to the path integral measure. As wc will 
ultimately see, the path integral approach to Ashtekar variables provides an 
infinite set of semiclassical states to choose from. A subset of these semiclas- 
sical states will correspond to the quantum states determined by the opera- 
tor canonical approach. Therefore, the requirement that the state satisfy the 
semiclassical-quantum correspondence will be equivalent to the requirement 
that the path integral and canonical approaches should be equivalent. 

It has not exhaustively been established the well-definedness of the path 
integral in metric variables as a means to corroborate the canonical approach 
to quantization of gravity. However, in ([11]) correspondence between ^path 
and ^can was shown at least at the linearized level of the constraints ap- 
plied to a linearization of the wavefunction using techniques similar to path 
integral quantization of Yang-Mills theory. One of the aims of this paper 
is to show that the path integral for quantum gravity is purely well-defined 
and convergent, as evidenced by the existence of the pure Kodama state. 




*[^] = <A(St)|A(So)) 

= j DAphysDAgauge^^V[iSAsh[A]] 

DApfiyg 



(10) 



where we have used (modulo metric signatures) 




(11) 
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2 Some more on the kinematic constraints 



This section is intended to present the kinematic constraints from another 
perspective. The reader can skip to the next section in order to avoid the 
details. The main purpose is to acquire as many tools as possible for our 
QGRA quantization toolkit. Three quarters of this section is intended for 
preparation for part II and III of this series. 

The the quantum physical states |*) = l^'can) are the set of states which 
lie in the kernel of the quantum constraints, whereas the semiclassical states 
I^VKfcft) lie in the kernel of their classical counterparts. 

For the kinematic constraints there is insufficient criteria to distinguish 
between these states, due to the fact that they are linear in conjugate mo- 
menta of the dynamical coordinate variables. However, the kinematic con- 
straints should be considered more fundamentally than consistency condi- 
tions on the model. Rather, they interact in a special way amongst each 
other and are intimately intertwined with the dynamics that uniquely fixes 
a state. Onr main theme of this series for each successive work, in addition 
to results, is to incrementally fill up our toolkit. 

We will explicitly solve these constraints for gravity coupled to matter, 
using the basis [4] for the CDJ matrix. The quantum Gauss' law constraint 
reads 



D, 



+ Qa] ^can = + Qa) ^can (12) 



where Qa is the eigenvalue of the matter part of the constraint on the state 
for matter fields (p"^ transforming in some representation Ta of SU(2), where 
the index A ranges from 1 to N, given by 



{G a) matter'^ can — (f>'^(Ta)BJ^^can — (4>^{Ta)B ^^B^^can — Qa^can 

^(13) 

and the Ansatz ^^an = has been used. The CDJ Ansatz a* = 

{^Ican/^AD = ^aeBl leads to the condition [4] that 

Xe^ae + Qa = (14) 

where the twisted vector field Xg is given by 

X,^ae = -^,^ae + {Cay'^fg. (15) 
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the SU{2)- 3-vector = te{x) can be thought of as a set of coordinates 
defining three directions in the tangent space (angles) of a nonorthogonal 
coordinate system defined by the vector field and corresponding covector 
one-forms 

^^ = 1^ = ^^^' dt^ = d^'(^~"yj (16) 
representing the rate of change of a function in any of the three internal 
SU(2) directions of the magnetic field. In this system we have (^dt^\di,^ = 6^ 
The new set of coordinates is related to the old by the relation 

f{x)=fiO)+ r dy\B-% (17) 

^0 

As long as is nondegenerate, which we expect in the presence of matter 
fields ([?],[4]), the map between the coordinate systems should be bijective. 
The integrability condition for these vector fields is given by the theorem of 
Probenius, which is the condition that there exist functions h^aj,{x) such that 

[^'|6]/(^) = ^a.(^)^/(^)Va,6 (18) 
for all smooth f{x). In this case the 'structure functions' are given by 

again which are globally well-defined due to the invertibility of the magnetic 
field. When the same criteria is applied to the 'twisted' vector fields, X^, 
twisted by a new kind of connection C = BA corresponding to a covari- 
ant derivative with respect to the new (nonorthogonal) coordinate system, 
taking values in a tensor representation of the group. 

{CaVa = BlA\{Uf5ge + febgSaf), (20) 

we obtain, in operator form 

[Xa, Xt,] = hl,Xe - K,Ce + Fab (21) 

where Fab is the curvature for the new connection Ce, a curvature of curva- 
tures, given by 

So it appears that in the presence of inhomogeneous terms, including 
matter, the ability of the new coordinate system to cover the spatial manifold 
S may be in jeopardy due to this obstruction. However, application of the 
theorem of Probenius enables one to transform into a new coordinate system 
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which does manifestly span S, even in the presence of matter fields, by 
finding an appropriate linear combination X = a^Xf, for suitably chosen 
coefficients a^. We will not attempt to illustrate this in this work, but 
will rather find it more convenient to adopt a fibre bundle interpretation 
on a base space G = {B~^)T, with fibres comprising the spin and spin 
2 elements of the CDJ matrix in the presence of a model-specific matter 
source, determined by the spin 1 elements and the charge. 

Recalling the notation from [4] ^ab = Dab + V'afc + ^ab where 

/ ^'13 

*a6 = *22 *23 

\ *3i *32 *33 

and the diagonal, off-diagonal and antisymmetric parts are given, respec- 
tively, by 

X00\ / U W 

Dab = \ Y \ ; ^l;ab = \ U V \ ; 

oozy \w V 

Using the notation, for integral curves along the direction t = f^, in analogy 
to path ordering of a parallel propagator, we define the 'symmetric contri- 
bution' D which acts on the symmetric elements of the CDJ matrix, given 

by 

DfcF^t'')^[U-'l-U]F, (23) 
where U is the parallel propagator along a path in this SU{2) space, given 

by 



U = U{t,0) = l- f dhCaibc){tl) + [ dh r dt2Ca(bc)itl)Ca{bc){t2) 

Jo Jo Jo 

rt rh i't2 
- dti dt2 / dt3Ca[bc){tl)Ca{bc){t2)Ca{bc)it3) + ■■■ 

Jo Jo Jo 

dt'C,^bc){t')]] (24) 




P 



exp 



and the antisymmetric contribution D which acts on the antisymmetric 
elements of the CDJ matrix, given by 

;^.abc Q 

De F{n^[U-'—U]F, (25) 
where the antisymmetric counterpart of the parallel propagator is given by 
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U = U{t,0) = l- dtlC„[i,,](ii)+ / dh dt2C^^hc]{tl)Calbc]it2) + ... 

Jo Jo Jo 

- dti dt2 / dt3Ca[bc]{h)Ca[hc]{t2)Ca[bc]{'t3) + ... 

Jo Jo Jo 

dt'Ca[bc]{t')\](26) 



= p 



with no summation over a, b, c, e, and with the definitions Ca{bc) = Ca6c + 
Cacb, and C„j(,(,] = Cabc ~ Cacb- Gauss' law constraint can be written in the 
form 



f)112 
^2 


Cl(23) 


] 


/ 1pl2 


f)212 




C2(13) 


1 V'23 


^3(12) 


n323 
-^2 


L>f3 J 


\ ■031 






^22 
^2 





n222 
-^2 







£>f3 



where Q' is made up of the matter charge and the antisymmetric part of the 
CDJ matrix, both determined by the matter fields, given by 





' D, C.23] 
^212 ^223 



,113 



Di 

V C'3[12] D2 



D3 (^2 [13]) 
A.323 A,313 





fA,2\ 


1 


A23 




U31 / 



Observe that it is in general nontrivial to solve for the shear components 
■012) V'23 and 013 in terms of the anisotropic components tpu, 022 and ^33 
due to the noncommutativity of the elements of the nondiagonal matrix of 
differential operators acting on the former. However, it is easier to solve for 
the latter interms of the former. 

The differential operators D and D have a useful interpretation. Sup- 
pose that a field F satisfies the following inhomogeneous first-order linear 
differential equation with source q{s) distributed along a curve parametrized 
by s 



d 



DF=[-+C]F{t)=q{t) 



(27) 



Then (27) can be viewed as an evolution equation in t space, which can 
readily be inverted to give 



011 

^^22 
0'33 



F{t, r, s) = U-\t, to)F{0, r, s) + U{t, to) f dt' U'^t', to)q{t', r, s) (28) 

Jto 
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(28) has a physically intuitive interpretation with regard to interaction 
among the CDJ matrix elements. The value of the 'fibre' F at any point 
along a curve t is determined by two contributions, (i) First, its value at 
the origin of the curve t = is parallel propagated to the spatial point in 
question {x{t), y{t),z{t)). This has a direct analogy in the Dyson formula in 
field theory, in which the parallet propagator plays the role of the interaction 
Hamiltonian. (ii) Secondly, there is a contribTition due to an inhomogeneous 
term which depends upon all points < t' < t. If this procedure can be 
carried out for three linearly independent directions t = {t^,t^,t^), then the 
value of F{t) can be found everywhere in E. 

This makes it possible to solve for the diagonal (anisotropy) components 
of the CDJ matrix in terms of the matter phase space variables, which are 
model-specific, and the symmetric off-diagonal (shear) components. The 
rotational (antisymmetric) components of the CDJ matrix are fixed by the 
matter contribution to the diffeomorphism constraint 



. ^^^^ 

where the matter contribution {Hijmatter = Hi is the eigenvalue of its re- 
spective operator on the state 



{Hi)matter'ifcan = Di^^ can = {Dicf)^-^^ ^ can = Hi^can, (30) 

which by the CDJ Anstaz requires that 



Aab = ^abdi>d = \H\ ^B^Hieabd- (31) 

Note that when there is no matter present, that the antisymmetric part of 
the CDJ matrix is identically zero. This is also the case for minisuperspace 
models not containing fermions. The relationship among the elements of the 
CDJ matrix imposed by the kinematic constraints, in general, holds for the 
state ^'can both at the classical and at the quantum level regardless of the 
model. It is the type and representation of the matter fields, if present, which 
distinguish one model from another. It is useful to visualize the diagonal 
(anisotropic) elements of the CDJ matrix as comprising a three-vector X = 
(X, Y, Z) and the off-diagonal symmetric (shear) elements as comprising 
another three- vector U = (U,V,W), with the antisymmetric (rotational) 
elements comprising a third 3- vector H = (u,v,w). The relationship among 
these vectors is linear and appears, in matrix form, as 



8 















X 
Y 

z 











n333 



f)112 
212 



\ ^3(12) 



Cl(23) 
n323 



A"' 



2(13) 









n222 
-^2 








n222 








^f33 




7f • (Ti,^) 

7f • (r3<^) 



,113 



i^i 



where the subscript denotes the value of the quantity on a 2-dimensional 
surface orthogonal to the integral curves of the vecor field comprising the 
original differential equation. Note, below that {t\^T2-,tj) need not form 
an orthogonal system of coordinates, but it should be possible to construct 
an orthogonal system from them using either Probenius' theorem, or the 
Gram-Schmidt orthogonalization procedure. 

^(0,Ti,r2) 
^(Ti,0,r2) 
Q , ^(0,0, ra) / 

Note that the first term on the right hand side, which involves the 
anisotropic components at the origin, can be expressed entirely in terms of 
the shear components and matter fields at the origin by solving the Hamil- 
tonian constraint. Hence 










222 





(i2, 




to))- 








( X[U{to)] 

y[u{to)] 

\ Z[U{to)] 



+... 



Making the identifications 



9 



D 



sym 













n222 
^2 


») 








J 


r)ii2 

^212 
^3(12) 


Cl(23) 
n223 

f)323 


C2(13) 
^313 



D. 



anti 



j^.112 

D2 
^212 

^1 

\ C'3[12] 



3 
323 



;:^113 



the final equation can be written, schematically, 



s s 

The interpretation of the canonically determined dynamics for the CDJ 
matrix can now be seen as the following neat picture. In order to determine 
the gravitational sector of the generalized Kodama state 



^ qrav CXp 



/ ^abBiAfl = exp f / {^abF" A - Ag(G„ - Q„))l (33) 



we need to find the CDJ matrix everywhere in M in order to integrate 
against the a 'topological four-form tensor' 

pa^pb^ ^^^^^p^^pb^dxi" A dx"" A dxP A dx" (34) 

which serves as a test function. For a given matter model the charge 
Q = (/>• (Tvr) and the momentum H = (f)-[DTT) are specified for all along with 
the contribution of the matter fields to the quantum Hamiltonian constraint 
0.1,0.2) (corresponding to the orders of singularity. This combination 
of inputs produces an ouput comprising the anisotropy X and the shear U 
components everywhere in space, and can be pictured as follows. 

Map 3-dimcnsional space S into an internal manifold G with coordi- 
nates {ti,t2,t^) by the action of a nondegenerate 'driebein' B*. Each point 
f of G consists of three SU{2) 3-vectors Xa = X = {X,Y,Z), Ua = 
U = (U,V,W) and cja = UJ = Co{Q,H). The pure Kodama state con- 
sists everywhere of two zero vectors Uxod = ^Kod = and the fixed vector 
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Xxod = -6(^GA)-i(l,l,l). These vectors remain globally decoupled ev- 
erywhere. For the generalized Kodama state ^' two vectors are specified: 
the 3- vector to' = u![Q,H], determined everywhere for by the given model 
(by compressing the matter fields into the gravitational arena [4]) and the 
shear Uq on three 2-dimensional surfaces (the {ti,t2), (i2,i3) and {ts,ti) 
planes, which is freely specifiable. 

The kinematics then takes over to give the anisotropy ^(i) everywhere 
in G, from which the shear everywhere U{i) = U[X{i)] automatically fol- 
lows via the dynamics. The coupling due to matter can then be seen as 
a transformation between two coordinate systems U U{X) (analogous 
to the transformation from Cartesian to spherical coordinates) induced by 
the freely specifiable quantities u and C/q. Each global configuration of 
{X, U{X)) corresponds to a generalized Kodama state. 

Since this relationship and interpretation holds at the semiclassical and 
at the quantum levels we must have that 1^*) = l^'vKfcb) with regard to 
the canonical relationship imposed by the kinematic constraints. Thus as 
the kinematic level there are three degrees of freedom, residing in the off- 
diagonal shear components of the CDJ matrix, that exist in the semiclassical 
state "if can = "^Wkb and in the quantum state "^can = and the states are 
thus far identical. The semiclassical-quantum correspondence is intact. 

Let us now establish the semiclassica-quantum correspondence via path 
integral for the simplest model, namely that of pure gravity with A term. 



3 The pure Kodama state: Equivalence of the canon- 
ical versrus path integral approach 

We would like to have an algorithm for determining the ground state wave- 
function corresponding to general relativity in Ashtekar variables coupled 
to the most general fields. First it is important to have an understanding of 
the pure gravity case with cosmological constant. Following suit as in the 
Hartle-Hawking prescription one can write 

^path = {A\A') = J DAexp{iS[A]) (35) 
where the action, by the 3-1-1 decomposition, is given by 



S[A]= [ (aiA^-e'^Ga-N'Hi-NH) = [ (alA^)-G{e)-Hi{N')-H{N) 

(36) 

Since the 3-1-1 decomposition is equivalent to transitioning from the tan- 
gent bundle to the cotangent bundle involves the phase space variables 
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{Af,A'l) {A1,al) the path integral should, strictly speaking, be per- 
formed over all phase space variables. Hence 



^path 



I 



DADaexp{iS[A, a]). 



(37) 



In a usual treatment the path integral of a constrained system with 
constraints <I>^ one would include a gauge fixing condition with a factor- 
ization over gauge equivalent orbits and correpsonding contribution to the 
path integral measure 



We shall, in our approach to path integral quantization, introduce the fol- 
lowing ideas: 

(i) We will not factor out any infinities due to redundant integration over 
gauge orbits, since these infinities will be common to all states and should 
therefore cancel out in the relative probabilities and in the computation of 
observables. 

(ii) As a consequence of (i) we will not need to implement any gauge-fixing 
procedures, and therefore needn't worry about the complications due to 
naive inequivalence among different gauge-fixing conditions. 

(iii) One possible interpretation of the path integral measure, if there is to 
be any equivalence between the path integral and canonical approaches to 
quantization, will be such as to enforce this equivalence. Since the differ- 
ence between '^can and '^path is, in our interpretation, a direct consequence 
of either the breaking or of the skewing of the SQC, which is due to the 
Hamiltonian constraint (a phenomenon unique to quantum gravity), one 
possibility (albeit artificial) is to define a measure 



dtx{a,A)^dtx{a,A)([[d{Ha{x))) 'llS{{^can)-'H^can) (39) 



which forces ^can = ^path providing the integral over the lapse density is 
performed prior to applying the measure. We will show how this measure 
can be usely to force by hand the equivalence of the canonical to path 
integral approaches, but in later work accomplish the same effect without 
it. But for now, let us focus upon the path integration over the phase space 
'coordinates' DA. 

When one decomposes the coordinate part of the path integral measure 
in (37) into a measure over the physical and the gauge degrees of freedom 



d(,{a, A) dtx{a, A) [] (5(x"(x))Det ([$^3, x"]+) • 



(38) 



X 



X 



X 




x,a,b 



(40) 
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one realises that the group averaging proccdTire [6] can only be performed 
with respect to the gauge group SU (2). Since the total phase space consists 
of the Ashtekar conection and the densitized triad a\ and excludes the 
auxilliary metric variables, qq^ and \/h, there is no inherent mechanism in 
(37) for implementing the corresponding Hamiltonian and diffeomorphism 
constraints via path integral. 

One way to implement the constraints would be to incorporate these 
quantities into the definition of the phase space variables and path integrate 
over them accordingly. This is the direct analog of integrating over 4-metrics 
instead of 3-metrics in the Hartle-Hawking prescription. The space of 3- 
metrics hij plays the role of and the qq^ parts can be viewed as structures 
external to the physical phase space which must be integrated over in order 
to reduce the state from an unconstrained form to a state that describes 
quantum gravity. In other words, to study quantum general relativity rather 
than simply SU(2) non-abelian gauge theory we must embed the phase space 
of the former within the latter by enlarging it. 

The configuration space of general relativity in metric variables contains 
6 — 3—1 = 2 physical degrees of freedom per point: 6 components of 
the 3-metric hij minus 3 diffeomorphism constraints minus one Hamiltonian 
constraint. The configuration space in Ashtekar variables contains 9—3—3— 
1 = 2:9 components of the connection A^ minus 3 vector constraints minus 
three Gauss' law constraints minus one Hamiltonian constraint (we have for 
the time being put aside the reality conditions). The number of degrees of 
freedom is the same for both but if we were to enlarge the respective phase 
spaces we would increase them by 4 for metric relativity (including 500 = 
and QQi = Ni). 

For the case of Ashtekar variables we would need to include and 
not N, but rather the combination TV = N/s/h. This is an increase of 4 
if we regard the lapse and the lapse density as independent variables, but 
5 if we do not. The former interpretation appears more reasonable, since 
we by definition regard any quantities involving the spacetime metric as an 
external structure with respect to the Ashtekar variables. However, although 
the shift vector iVj is external to the 3-metric hij as it is to the Ashtekar 
connection A1, it would not be obvious a-priori to conclude that the lapse's 
externality to the 3-metric is the same as the lapse density's externality to 
the Ashtekar connection. One can say that they must somehow be equivalent 
if the two theories are to in a sense correspond to each other. But as we 
can already see, the Ashtekar variables 'know' something about the metric 
general relativity irrespective of the quantization scheme. We nevertheless 
take the Lagrange multipliers N_(x) and N'''[x) as external structures in the 
reduced theory. 

To implement the Hamiltonian and diffeomorphism constraints we must 
do so by hand. The most seemingly natural way is to act on (37) with the 
unit operator in the form 
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f DNDN' , , 

/ = T ^- (41) 

This has the effect of implementing the constraint while introducing a con- 
tribution to the normalization factor of the wavefunction. Since we know 
a-priori the right answer, namely the Kodama state we can place restrictions 
upon the form of the normalization factor. Thus we have 

*[A]= (^J DNDN'Y^ J DNexp[iH{N)] J DAfDA^exp[iG{e)] 

X j DN'exp[iH\Ni)]exp(i J a\M^. (42) 
where we have introduced the following notation for smearing a function, 

F'iGj) = Y.[ d'xF{x)G{x) (43) 
. Jm 

Alternatively, if we adopt the enlarged phase space interpretation then we 
can discard the prefactor entirely. The first two integrations transform the 
unconstrained wavefunction 

^o[A] =exp\i [ aiA'}] (44) 

into a gauge- invariant, diffeomorphism- invariant functional via a group av- 
eraging procedure [ ] 

^inv[A]= J DNWe''^[N\9'',A] (45) 

where 



*[iV^r,^] = exp[iHi{N')]exp[iGa{0'')^o[A]exp[-iGaexp[-iHi{N')]. 

(46) 

Here 'i'inv [A] is a diffeomorphism and gauge invariant state derived from the 
unconstrained state and is otherwise arbitrary- this is the maximum that 
can be done at the the kinematical level. Also note that (46) expresses the 
tranformation as an operator relation. This is possible because quantum 
operators corresponding to the diffeomorphism Hi and the gauge Ga con- 
straints are linear in functional derivatives and therefore can be expressed as 
eigenvalues directly in terms of a phase without quantum corrections. This 
property is not spoiled due to ordering ambiguities due to the fact that the 
kinematical algebra of constraints Diff * SU (2) closes, as a lie algebra, on 
itself. By the Baker-Campbell-Hausdorf formula 
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(47) 

one ends up with a linear combination of diffeomorphisms and gauge trans- 
formations in the exponent which is still first-order in the momenta. As 
the Hamiltonian constraint commutes with the Gauss' law constraint but 
not with the diffeomorphism constraint (except for spatially uniform lapse 
density functions N_) the question arises as to whether the sequence of per- 
forming its path integration is important. 

As in any theory given possible ambiguities due to ordering one chooses 
a particular operator ordering and then quantizes the theory with respect to 
the chosen ordering. Then one can simply transform among states correp- 
sonding to different ordcrings. Here, we will isolate the kinematic portion 
of the phase space first by performing the path integrals as above, and then 
apply the Hamiltonian portion last in order to derive the physical state of 
the model. 

That the starting wavefunction ^'o[^] is not a-priori gauge or diffeo- 
morphsim invariant can be seen by performing a transformation 



*o[^] = ^[A, a] = exp[if{a)]exp(^ J a* if ) exp[-zf (a)] 

(48) 



expi 



M 



where T{a) represents either a SU{2) (g) Diff 'gauge' transformation la- 
beled by a = {9"',N^) as determined by the BCH formula. We are re- 
stricting, for simplicity, to 'gauge' transformations connected to the iden- 
tity. In any event, the result is to produce a wave- functional that satisfies 
the Diff~kSU{2) constraint, since the path integral can also be viewed (by 
'disentangling' the parameters) as 

= \{5{Ga{x))6{Hi{x))^o[A]. (49) 

x,i,a 

which is the statement that the state has support only on configurations 
for which the kinematic constraints are satisfied. The projector approach to 
path integration of constrained systems is well explained in [6]. A special 
case is that the kinematic constraints are identically satisfied when the self- 
duality condition (cr* = kB^^ for some numerical constant n) holds. Imposi- 
tion of the Hamiltonian constraint fixes the value of this numerical constant. 
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To implement the Hamiltonian constraint via path integral we can attempt 
to define a state evolved in time relative to '^inv via 

^{A,N) = exp[ii?(iV)]'I'i„,[^] (50) 

and then integrate over all N_. The difference from the previous procedure 
is that the Hamiltonian constraint is cubic in momenta, which introduces 
an operator ordering ambiguity, as well as potential divergences. Since the 
path integral deals with c-number quantities and not quantum operators, it 
would not be capable of dealing with the Hamiltonian constraint other than 
at the classical level. 

One aspect of the ^can = ^path equivalence is that consistency between 
the two must be maintained at all stages for the particular state in question. 
The CDJ Ansatz clearly picks out a range of states canonically. One would 
expect the path integral to reflect this and can taylor it accordingly. For 
an ordering with momenta to the left of the coordinates, one has in the 
canonical approach under the Ansatz ^can = e^''"" 



*can = (e.jfcT7i^^a) *can = (51) 

By the CDJ Ansatz a* = (blcanl 5A1^ = "^aeBl we find that the antisymmet- 
ric part of '^ae must be identically zero. Hence the CDJ matrix is symmetric, 
both at the classical and at the quantum levels. This is the case due to the 
linearity of this particular constraint in momenta 

For the quantum Hamiltonian constraint, for the chosen operator order- 
ing of momenta to the left of the coordinates, this yields 



B'l + hG- 



6 \5Al. 



*can 



Note that the quantum Hamiltonian constraint is a scalar constraint, and 
therefore in the general case can impose only one condition upon the phase 
space. A nontrivial solution can be found such that 

Bl + hG^(^-^)=Q^k,c, (53) 

which is tantamount to converting the cubic Hamiltonian constraint into a 
linear constraint, upon the same footing as the kinematic constraints. This is 
as well appropriate for the path integral, since the aforementioned infinities 
and ordering ambiguities are no longer an issue. By the CDJ Ansatz one 
has 
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(*ce + 6{hGA)-^6ce)B^ = (54) 

One obvious solution to (54) is that ^ce = —6{hGA)~^Sce [2] which, is 
consistent with the Gauss' law constraint and leads directly to the pure 
Kodama state 

^Kod0^e-^(~'^^r'ics[A]_ (55) 

This corresponds to nonflat Ashtekar connections. From the perspective that 
the Hamiltonian is a scalar constraint and Gauss' law is a SU{2)^ vector 
constraint there must remain, in the most general circumstance, two degrees 
of freedom in the state. Noting that (54) is in reality a matrix acting on 
the SU{2)^ index of the magnetic field B^, viewed as a collection of three 
3-vectors labelled by the spatial index k, the condition that the equation 
hold for nontrivial is weaker than that leading to ^Kod- One can also 
have 

det (^*ce + 6{hGA)-^dce) = (56) 
which leads to the requirement that 

det* + ISA" Var^' + 216A-2tr^' + 216A~-'^ = 0. (57) 

(57) is one equation in five unknowns. = ^6{hGA)^^Sc.c is a special 

solution, leading to the pure Kodama state, but it is clear that there are 
others as well. Hopefully the path integral is 'smart' enough to pick out all 
solutions provided by the canonical approach. We will demonstrate this in 
the next section. 

Note that since ^ab is symmetric it can be diagonalized by orthogonal 
transformation 

^ae = OabDbcOl (58) 

where Dab = Diag{X,Y, Z) is a diagonal matrix and Oab = O{9i,62,0-i) is 
a 3 by 3 orthogonal matrix parametrized by three angles {9i, 02,03)- Note 
that all quantities appearing in (57) are independent of Oab since 

det* = (detOfdetD = dctD = XYZ 
tr* = Ol^OacDbc = tiDbc 
tr*2 = Oab{DD)bfO% = trZ?2 

(59) 

Hence Var^ = VarD. This leads to the condition 
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XYZ + 36A'^{XY + YZ + ZX) + 216A~^(X + Y + Z) + 216A~^ 

V 216 + 36A(X + y) + y ^ ' 

which determines the diagonal (anisotropic) elements of the CDJ matrix 

in terms of the two freely specifiable parameters X and Y . This leaves 
the Gauss' law constraint, three conditions upon the off-diagonal symmetric 
(shear) components, which should be specifiable completely in terms of the 
diagonal ones. 

The CDJ Ansatz can be viewed as the imposition, by hand, of a con- 
straint upon the path integral which collapses the phase space path integral 
into a configuration space path integral. Returning to the Hamiltonian con- 
straint portion of the path integral we have, in an abuse of notation, 

*//am[A?] = / Dm!{A,a,N) = ]\5{H{x))^inv[A,a]. (61) 

At the level of (61) all constraints have been imposed subject to the CDJ 
Ansatz, which produces a two-parameter space of solutions for the chosen 

operator ordering, for which the Kodama state is a specialized solution from 
this set. Since the semiclassical-quantum correpsondence is maintained we 
have 



{a,A){\{5{H,l{x)))-'\{5{{^,anr'H^can) = = 



(62) 



and the path integral becomes, taking into account the substitution of the 
CDJ Ansatz, 

"^path = J DADa n ^i'^iix) - ^aeBl{x))^o[A] (63) 

x,i,a 

upon path integration over the 'momenta' the Ansatz is substituted, 
along with the solution imposed by the the canonical constraints, into the 
remainder of the path integral yielding [1],[2] 

%ath = ^pathiX, Y) j L>^exp [j^ ^F'> + AlDi{ai)) 

= J DAe^p [J^ [^ab{X, Y)F'' A F^] (64) 

To enforce the Gauss' law constraint we shall illustrate via the method 
introduced in [1], for simplicity. The Gauss' law constraint can be written, 

in the absence of matter, 
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Dial = di 



0, 



(65) 



where in (65) we have used the path-ordered parallel propagator along a 
path parametrized by s, given by 



J'y J sr\ 



dsxWaAf yU = Pj~<^ 



(66) 



It is clear that the pure Kodama state is a particular solution to the 
Gauss' law constraint since Qa = when there are no matter fields present, 
since 



-^{hGK)Xe5ae = -Q{hGK) (^ — Sae + fahcA\^ce + /e6c4*ac) 

= -Q{hGK)-\fab^A%e + febcA%c) = (67) 

due to antisymmetry of the structure constants. However, there are more 
general solutions, which correspond to more general states. Applying the 
abelian Poincare Lemma and the CD J Ansatz to (65), 

(PeA^)^a^ = (PeA^) Vfeei?:) = e'^'djVka. (68) 

for some arbitrary SU{2)^- valued 3-vector Vka- Solving for the CDJ matrix 
we have 

*6e = {U-^rb{djVka)e^''\B-^)ie = He (69) 

At this point Vbe is an arbitrary S'J7(2)_-valued 2-index tensor. The CDJ 
can be expressed in the in the spin-2 representation in terms of the rotation 
angles 

^ae{0) = Oab{e)D,,Ol{e) = Rt{6)Djy = (^"1)^ (^1 , (70) 

where Dfg is a diagonal matrix. Note that since ^'^e is symmetric, three de- 
grees of freedom can be eliminated from Vae by symmetrizing on the indices 

= {U-%Vae). (71) 

The diagonal elements of the CDJ matrix can then be expressed, with the 
rotational degrees of freedom 9 isolated, 

Df, = {R-\e))%{U-'V),,. (72) 

(72) can be compared with the form derived in sec(2), in terms of degrees 
of freedom 
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(73) 



In this basis Xq is the anisotropy vector at the origin, which corresponds 
to the diagonal matrix elements, which fixes three components of Vae- The 
remaining three components of the latter must be chosen to vanish in order 
to balance the degrees of freedom. The shear vector U corresponds to the 
degrees of freedom in the rotation angles 9. 

Let us now see how Gauss' law interacts with the Hamiltonian constraint 
to determine the final CDJ matrix elements. U can be considered chosen at 
all space and Xq is arbitrary. We will discretize space, labeling the points 
by subscripts, to illustrate. Since the Hamiltonian constraint is independent 
of 9 it is the former which determines the latter at each stage of the process. 

(i) (Xq, Iq)) {U{j, Vq, Wq) are freely specified, (Xq, Yq) being substituted into 
the Hamiltonian constraint, which fixes Zq. There are five degrees of freedom 
at the origin. 

(ii) Using (73) and Ui as an input, Xi is found explicitly as a functional of 
(Xq, lb) Uq, Vq, Wq) and (C/i, Vi, Wi) producing the vectors 

Xi = Xi{Xq,Yo,Uo,Vq,Wq,Ui,Vi,Wi) 
Yi = yi(Xo, Fo, f/o, ^0, VFo, f/i, Fi, Wi) 
Zi = Zi(Xo, Fo, Uq, Fo, Wq, Ui,Vi, Wi). (74) 

(iii) Consistency must of course be checked with the Hamiltonian constraint. 
Sustitution of (74) into the Hamiltonian constraint reveals that not all com- 
ponents of C/ = {Ui,Vi,Wi) are independent. This limits one of the com- 
ponents, say W\, to be a functional of all the variables thus far, including 
{Ui,Vi). Thus there are only two degrees of freedom accepted from U. 

(iv) The new vectors are substitued into (73) and the process repeats until 
all of space is covered. 

To summarize, the CDJ matrix elements is freely specified at the origin 
except for the antisymmetric components, which arc zero due to the absence 
of matter. These components remain zero throughout the 'evolution' (here, 
evolution with respect to space, not time) . Subsequently, two components of 
the shear are freely specifiable throughout S, from which the third compo- 
nent and all anisotropy components are uniquely determined. This should 
yield a two-functional-parameter family of generalized pure Kodama states. 

We have found the most general state when there is no matter present, 
even more general than the pure Kodama state, due to the two aforemen- 
tioned degrees of freedom. When there is matter in the model, the matter 
charge and diffeomorphism momentum are injected at each stage in the 
process, which can be considered a nonlocal effect. 

The pure Kodama state is given by the 'initial' conditions 
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{X, Y, Z, 01 , 02, e^)Kod = -6{hGA)-^ (1, 1, 1, 0, 0, 0) (75) 

for pure gravity in the absence of matter, whereupon a vanishing shear input 
at each stage maintains its CDJ matrix elements globally decoupled. This 
would not be the case for the general matter model, which is labeled by 
a two-parameter family of functions corresponding to the shear elements 
U{x),V{x). 

Note that for the pure Kodama state we have "^ab = — 6(/iGA)~^5(jb and 
(64) becomes 



path 



= f DAexp\-6{hGAy^ [ trF A f] (76) 



and something interesting happens. The path integrand becomes purely 
topological. By Stokes' theorem, 

/ tiFAF= [ tr{AdA + lAAAAA)=Ics[A{^)]. (77) 
Jm JdM ^ 

The path integrand of (76) collapses from a volume into a boundary integral, 
which becomes immune to path integration within the interior of M and can 
be factored out. 

'^path= [ DAe^p\-6{hGA)-^ [ tiF A f] = [ DA^Kod[A{^)] 

= (Fo/A)exp[-6(nGA)-i/c5[^(5])]] (78) 

The Volji factor, the volume of the space of Ashtekar connections, is 
formally infinite. But we expect this infinity to cancel among ratios of 
states and for observables. The result is a topological state, defined on 
the boundaries of M, just as the Hartle Hawking state is parametrized by 
the boundary values of the 3-metric hij. The path integral resulting in the 
Kodama state is manifestly convergent. We expect to be able to perform the 
analogous procedure for a general model when there is coupling to matter 
fields involved. 

We will now show a sense in which the states ^Kod{U, V) are orthogonal, 
for this and in preparation of future work in the series. Note also that these 
states are quantum states, even though we solved only the classical part of 
the Hamiltonian constraint. This is the case becasue the state, as ^A-od[^]) 
is determined for an operator ordering for which gi = §2 = a-priori. 
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4 Generalized Kodama wavefunctions as a com- 
plete basis of states 



We have demonstrated in the previous section that there exists a two- 
parameter family of states satisfying the quantum constraints of pure grav- 
ity, in addition to the pure Kodama state formed by the Chern-Simons 
wavefunction. The first matter of interest is the sense, if any, that these 
states form an orthonormal basis. Secondly, assuming that there is such a 
basis of semiclassical/quantum states- which respect to what measure? The 
same questions will be raised of the more general models coupled to matter. 
We will demonstrate in this section the fundamental concept, which we hope 
partially addresses the issue of normalizability of the states. 

In order to establish in what sense an equivalence exists between the 
canonical and the path integral approaches to quantization and the normal- 
izability of the states that they produce, we must introduce a new repre- 
sentation for the Hamiltonian constraint. Wc will motivate the formalism, 
drawing upon simpler examples of increasing complexity, starting from the 
finite dimensional analog. 

Consider a polynomial F{x) which has A'^ zeroes, ri, r2,...riv- The poly- 
nomial can be written in the following factorized manner 

N 

F{x) = l[{x-rk). (79) 
k=l 

(79) is written based upon a finite, discrete spectrum of roots. This can be 
extended to transcendental functions with an infinite number of discretely 
spaced roots, for instance in the Euler expression for the sine function 

oo 

sinx = xYl{x^ - n^TT^). (80) 

n=l 

To gain an intuition for what is to follow, let us first extend this concept to 
the case of a continuous spectrum of roots. The Dirac delta function S{x), 
along with its various representations, can be viewed in a certain sense as 
a function with a continous spectrum of zeroes, namely the entire real line 
excluding possibly a finite number of points. One may imagine an extension 
of (79) for which the roots become closer and closer spaced together, and 
taking the continuous limit 

S^{x) = lime^o^^Y\.ix — Vk) such that||rfc — rfe_|_i| < e. (81) 
k 

or variations theoreof. 
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While we are discussing such pathological, albeit extremely useful, 'func- 
tions', let us generalize the concept of roots to functionals on infinite dimen- 
sional spaces. 

Bear in mind that what follows is meant, not so much to be a rigorous 
mathematical construction, but rather a qualitative argument to motivate 
the plausibility for how and why the the states satisfying the constraints can 
be in a sense be considered normalizable. 

Consider first the classical Hamiltonian constraint H{U,W,V) = Vx. 
Let W = W^'^\U, V) be a set of N solutions labeled by (7) for each U and 
V. The parametrization of the Hamiltonian constraint must reflect the fact 
that there is as well a two-parameter familty of solutions (or roots, if you 
like) U = C/(°^ and V = for each a = p. The classical Hamiltonian 
constraint function, which is explicitly a function of the three remaining 
CDJ matrix elements, can then be written as 'functional' factorization over 
each solution. 

Let us now discretize space, assigning the label n to each point. The 
Hamiltonian constraint, expressed in this language, reads at the classical 
level when 



n a 13 

<^(n(^(^-) - W^'^\U^'^\xn)M^\xn)))) (82) 



For brevity we we will perform the analysis for one spatial point and the 
compose the results in the end in for all points in the continuum limit. It 
makes no difference to our analysis whether a and /3 are discrete or contin- 
uous indices. Think of a and (3 as indices in functional space, each index 
labeling a different pair of functions for which the HAmiltonian constraint 
is satisfied Mx. In order to aid the visualization of the mechanisms at play, 
we shall explicitly maintain the indices throughout the derivation. Using 
the identity 



where the roots of F{x), we have for our discretized system of points 

and functions. 
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(83) 
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Ar(a)_l 



o. 71 



7 737^7 



-1 



,(84) 



where Ar("^), M^^^ and iC^''') are meant to take into account any possible 
multiplicity of 'roots'. In shorthand motation, 



Q!,/3,7 

5{V{Xn) - V^f'^XnmWiXn) - W^-'Hxn)) 

with a Jacobian contribution to the measure given by 



(85) 



f^a/3j{Xn) = iV(")(x„)M(/^)(x„)KW(x„) X J] {U^^H^n) ' U^^^^Xn))' 

71^0,72 7^/3,73^7 

{V(^\x„) - vM^x„)f''-\w(''\x„) - T^(-3)(,^))iv(^'-^ 

Again, a and /3 can be viewed as coordinatizing a 2-dimensional surface 
([/("), in the functional space over each point x, upon which various 
quantities of interest can be defined. 

We wish to evaluate a path integral for a functional ^' [U, V, W] upon 
which we wish to impose the Hamiltonian constraint. Hence, 



%ath = j DU j DV j L>Ae^^W*[C/,y,W^], 
where the path integration measure is given by 

DUDVDWDX = Y[dU{x)dV{x)dW{x)dX{x) 

X 

n 

= lirrin^oo JJ dU{xk)dV{xk)dW{xk)dX{xk) 



(87) 



(88) 



Note that *patfe is at this stage a single wavefunctional. Now observe what 
happens. First, we perform the path integration over the Lagrange multi- 
plier to obtain 
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%ath = / DUDVDWl[6{Ha{x))^{U, W, V). (89) 

Before proceeding with the path integral we must first distentangle the con- 
straint into a form Unear in the variables of path integration, including the 
appropriate Jacobian factor (the analog of (83) , except in the space of func- 
tions). We anticipate that imposition of the constraint will fix the function 
W for each set of functions U and V, which can be chosen arbitrarily. Thus 
we have, in disretized notation. 



%ath= Yl / DUDVH / dW{Xn)v^''^[U{Xn),V{Xn),W(Xn)]{f^aPyr\xn) 

S{W{Xn) - W(-''^[U{Xn),V{Xn)mU{Xn) - U^^'HxnmV (Xn) V^^\xn)) 

^U{Xn),V{Xn),W{Xn)],m 

where r]{xn) is the Jacobian of the distentanglement relation for W . Since 
the classical part of the Hamiltonian constraint is already linear in W the 
relation is given by 



^(7) = (;7"y.«)-i 

W [U[Xn),V[Xn)l)- A 216 + 36A(C/W + F(/3)) + A^C/W^/^) J-^ ^ 

Observe that the wavefunction, which was originally a single state, has 
now decomposed into a linear combination of states. There is one state 
for each pair of functions U{x) and V{x), consistent with the Hamiltonian 
constraint. We have taken into account all possible functions, and therefore 
all possible states. 

First we perform the W path integral, whereupon the three-dimensional 
functional manifold over each point collapses into a two-dimensional func- 
tional manifold. Now the path integral shifts to DU. This can readily be 
visualized by undiscretizing after the DW integral and the re-discretizing 
in preparation for the DU integral. Here, the notation ^'(2)) signifies the 
collapsed path integrand, now containing two functional degrees of freedom. 



a,/3,7 ' 



%ath = E / DVDUr,^^^l[U{x),V{x)]{fi^fSj{x))-'6{U{x) - C/(°)(x)) 

(5(F(x) (a;))^' (2) [[/(x),F(x)] 

iTTi^ ^ T. 

/(2) 



/ DV]JdUiXn)r]^^hU{Xn),V{XnWaP'y)-^S{U{Xn)-U("\xn)) 



a,/3,7 n 

S{V{Xn) - V'^P\xn))^^2)[UiXn),V{Xn)].i92) 
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The path integration over DU imposes the delta functional giving 

%ath = E / DVv^^][U^''\x),Vix)]{i,o,p^)-'[U(''\x),Vix)] 

S{V{x) - F(^) (93) 

The wavefunction in the path integrand has now collapsed to one functional 
degree of freedom. Again, we undiscretize and re-discretize to do the DV 
path integral, 

^path = En/ dViXn)vlll[U^''Hxn),V{Xn)]itiap-yr'[U^''Hxn),V{Xn)] 

5{V{Xn) - y(^^(Xn))*(i)[f/(")(Xn),y(Xn)] 

So the state has collapsed into a linear combination of the set of all functions 
on a reduced phase space consistent with the constraints. Furthermore, these 
functions are weighted by their respective Jacobians. To ease the notation, 
one can regard either a and P or U and V as dummy labels, with the 
understanding that each pair of indices labels a particular functional of two 
variables (functions) of the set. 

This procedure can in principle be carried on for more sets of variables. 
Now we must ask in what sense the states are normalizable. Let us say that 
the there were additional variables A in the path integral that needed to be 
done. A possible trap to fall into is to define the norm by 

<*a„ft|*a,,/3,> ^ / DA^^^,^^[A]^^.,f,.[A]. (96) 

(96) is incorrect because it is neither the correct set of variables nor correct 
Hilbert space with respect to which the orthonormality naturally would 
be assessed in the analogous system of ordinary quantum mechanics. To 
provide a more familiar analogy, imagine that a wavefunction ^ expanded 
in a complete basis of states 

^{x,y) = ^CnMx,y) (97) 

n 

where c„ are the weights. Performing (96) would be analogous to defining 
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(V-mlV-n) = dy U dx^m{x,y)^ * ( / dx^nix,y)) (98) 



where the product * would have to be suitably defined. The proper defini- 
tion, in accordance with the axioms of quantum mechanics, is given by 



with due regard for the line bundle structure over the {x, y) base space. The 
Hilbert space is the set of square-integrable functions on the the space of 
positions (x, y). In the case of the path integral the role of x was played by 
U{x) and V{x), variables which have already been integrated out, and the 
role of y played by the remaining variables, such as the connection A1{x). 
In order to assess the inner products and norms of the functions, one must 
integrate the probability density with respect to the variables in question: 
square then integrate, not the other way around. 

Applying this principle, we would have to go back to the corresponding 
stage upon path integrating W and isolate the individual states, including 
their measures and delta functions. It so happens that the states are peaked 
on functions consistent with the solution of the Hamltonian constraint. 

The proper state with respect to which orthonormality should be mea- 
sured is a functional of U{x) and V{X), labelled by a pair of functions 11°^ 
and F". So we have (suppressing the suffix on ^path) 



<*„,;3|[c/,F]) = 4c/ <5[F - F(^)]p[[/("),F(^)]*(2)[[/,y]. (100) 



Assuming a complete set of functional basis states, the the inner product of 
two states is given by 



[V - [u - [V - ^(2) [u, v]^^^)[u, v] (loi) 




(99) 




Tidying up the notation a bit, we have 
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DUDV5[u - [V - [u - [/(°')] [V - 

= JJ^[C/{aO _ u(-')] [V^^') - |p(a«,/3»)|'|M/(2)(«»,/?»)f . (102) 

In a nutshell, we have found that the inner product of two different 
generalized Kodama states, corresponding to any two different solutions to 
the constraints, is zero unless the two functions paramterizing the two states 
are identically equal. This it the orthogonal basis we seek: orthogonality of 
functionals with respect to arbitrary functions. To calculate the norm of a 
state we merely 'contract' the 'indices' and obtain 



(103) 



It may appear that (103) is pathological, due to the infinite product of delta 
functions of zero, two per spatial point. However, we expect these infinite 
numerical factors to cancel out in the ratios of probabilities amongs different 
states and also in the computation of observables. The functions of course 
must all have the same domain in order for the relative probabilities of 
solutions to the constraint in a superposition state to be well-defined. Still, 
( function rcgularization procedures can be used on these delta functions, 
although not necessary. 

This entire argument is based upon assessing the orthonormality with 
respect to the appropriate Hilbert space. There is of course still the possi- 
bility to carry out this computation in the remaining Tinintegratcd variables. 
The question is then to find and to to appropriately define with respect to 
what relationships orthogonality is being measured. 

Y[ J dU{x)AdV{x)AdW{x) j DN/"^^ dU AdV AdW5{Hd{U,W,V)) 

(104) 

which, when disentangled and account is taken of the state, leads to 



Y[dUdVdW {nap-yT^SiU - [/("))(5(F - V^^'^)5{W - T^^^^) 



exp 



M 



^ae{U,V,W,u)BiA-^+T(4) 



5^(At„;3M)-iexp[/ Uae{U^''\V^P^)BiA1 + T(.^)\ (105) 
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For brevity we have suppressed the dependence upon W and 7 above, but 
we can see that at the level of implementation of the quantum Hamiltonian 
constraint, the original path integral has split into a linear combination of 
(uncountably infinite) basis states. 



5 The problem of time in quantum gravity 

Let us revisit the no boundary prescription for determining quantum gravi- 
tational wavefunctions [8] . The no-boundary proposal relates the transition 
amplitude for the quantum gravitational wavefunction into its present state, 
to the path integral via 

{hij{^T),<P{^TMj{^o),'P'i^o))= [ DgDcPeMi [ LEnigA]) (106) 

J J M 

where the gravitational portion of the path integral is over all 4-metrics 

g^^ throughout the interior of a 4-manifold M with the given 3-metrics 
and matter fields (hij{YiT)-,4>{^T)) and (/lij(So), 0(So)) on the spatial 3- 
boundaries and Sq. The Lorentzian path integral, as written, is typically 
ill-defined. 

The usual prescription to find the ground state for a quantum field theory 
field theory is to Wick rotate the time interval t = T — —ir into 

the imaginary time axis in order to improve the convergence properties of 
the path integral. To isolate the ground state of the system one inserts a 
complete set of orthonormal eigenstates of the total Hamiltonian operator 

I = ^\n){n\. (107) 

n 

Thus 

{hij,<t)\h'ij,(t>') = J2 -^n[hijAWn[hij, 0']exp[-E„r] 

n 

= J DgDct)eM-S[g,(t>]), (108) 

where En is the n*'* energy eigenvalue of the quantum Hamiltonian operator. 
Then one takes the limit t ^ 00 and observes that it is only the n = term 
which survives, n = corresponds to the lowest energy level of the system 
and relates to the amplitude to make a transition from a time far in the 
past, to the present when measurements can be made on the system. 
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When this prescription is extended to quantum gravity one is faced with 
our interpretation of the problem of time. Succinctly, in order for the ground 
state to be selected, it must be incapable of having evolved into the current 
state, due to the condition that the state must identically satisfy the Hamil- 
tonian constraint. The requirement of equivalence between the canonical 
and path integral approaches for a reparametrization-invariant theory such 
as general relativity, to quantization system seems to be a root cause. Let 
^can and ^path be the wavefunctions determined, respectively, by the canon- 
ical and the path integral approaches to quantization. 

The quantum Hamiltonian constraint reads 

H\^can)=0. (109) 

Hence for compact manifolds or for manifolds with suitable fall-off conditions 
of the fields, it seems that there can be no evolution of the quantum state. 
Written in the language of transition amplitudes, 

{hij{J:T),H^T)\hij{^0),4>{^0)) 

= ll6{hij{ET) - /ii,(Eo))<5(</)(ST) - <^(So)). (110) 

X 

The no-boundary prescription circumvents this issue by restricting the 
initial state of the universe to the hypersurface t = —ioo, [8] for which it can 
be argued that the state does not exist. To entertain the possibility of the 
current state of the universe having evolved in finite time from something 
other than nothing, we shall focus on two spaceial categories of states: (i) 
There is a canonical state ^can^ in the full theory of quantum gravity which 
must satisfy the quantum version of the Hamiltonian constraint 

H\^can)=0. (Ill) 

For the sake of convenience let us assume we are now describing Ashtekar's 
gravity, although the same argument can essentially be applied as well to 
the metric representation to a certain extent) . (ii) But there is also a semi- 
classical state also of the full theory, "^wkb, that satisfies a condition usually 
of the form [2] 

H\^wkb) = {Ha + H,t}\^wkb) = {hG6^^\0)q, + h^G\6^^\0)fq2)\^wkh)- 

(112) 

where qi and q2 can be expressed in terms of the CDJ matrix elements in 
the form 



Qi = ^ijk^ 



abcG^h 



1'ae + ^^2 + 72 = (113) 



30 



+2hGA{5lC) - SlcD^ae'^bf + J^i = 0. (114) 

For this state, the semiclassical part of the Hamiltonian vanishes Hd = 0, 
but there remain quantum terms Hd- It is these quantum terms that can 
potentially provide a 'Hamiltonian' for the evolution of the semiclassical 
state. The perception of time evolution for this state is the classical 
limit of an unobservable quantum effect. So the first question that needs 
to be answered is, which version of the constraints is being satisfied and by 
which state? Since we make contact with the semiclassical limit of quantum 
mechanics, one may think that when we observe time evolution classically, 
we lose the ability to do so quantum mechanically, hence that it is impossible 
to experimentally verify any effects of quantum gravity, although we can 
verify the classical limit of such effects, such as inflation, etc. 

Part of the reason for this paradox may reside in the notion that one must 
choose the initial state, at — ■ioo, to be the only one with respect to which a 
measurement of the present can be made, in order to select a ground state 
consistent with the Feynman path integral. The no-boundary prescription 
would allow for excited states which when unsuppressed can tunnel into 
our present universe. Such states, evlauated at finite times, can interfere 
with a quantum mechanical measurement made at finite times. However, 
an alternative is suggested. 

One would like to be able to observe quantum gravitational effects exper- 
imentally, which requires a meaurement of two observables separated by a 
finite time separation r. The ability to correlate such measurements would 
be strong indication for the testability of the quantum theory of gravity 
below the Planck scale. 

First, a more convenient basis is neededed for measuring energy than 
the mode eigenstates |n). Let us consider instead a set of complete states 

quantum states [^can), which are solutions to the quantum Hamiltonian 
constraint. Wc will later have to define a measure with respect to which 
these states are complete, but assume for now that there exists a resolution 
of unity of the form 

^ = El*-n><*il|- (115) 

n 

The index n in (115), which can be discrete or continuous, enumerates these 
quantum states. Now insert these states into the quantum amplitude to 
transition from one state at time ti to a state at time tf for a finite time 
separation ti —t2 = r. The transition amplitude then reads 
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n 

n 

n 

J JM 

In the second line we have made use of the fact that the kinematic constraints 
have been identically satisfied, and that the Haniiltonian constraint has been 
satisfied as well, and we have used the smearing convention 

H{N) = [ (fxN{x)H{x). (117) 

JM 

Note that in this basis we have related a ground state, or a linear combina- 
tion of ground states, of quantum gravity to the path integral. Furthermore, 
these ground states occur at a time interval apart which should make possi- 
ble a measurement and correlation of observables. There is no interference 
from excited states because all states in this set are by definition of zero 
energy. 

It is a matter of convention whether one chooses a Euclidean or Lorcnzian 
signature for spacetime, but either way, if there is a quantum state satisfying 
the constraints, then it must necessarily be related to the path integral. 
Therefore the path integral and canonical approaches are still congruous. 

Now, let us assume that there exists a complete set of semiclassical states 
permitting a resolution of unity of the form 

n 

Again, the sense in which these semiclassical states from a complete set will 
be discussed future work. Observe what happens when we insert this set of 
states into the transition amplitude. 

E<^?(St,),0(S*,)|<l,,><*H|^f(S,J,</,(S,J> 

n 

= i: *£L[A"(St,), </>(S,,)]e^^^(^')e^«^^(^")*i;)LK(S*J, 0(S,J]e-^c.(iV) 

n 

n 

=? / DAD(l>eMi f LAsh[A,mm 

J JM 
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Here we have made a few assumptions, (i) The kinematic constraints have 
also identically been satisfied by the WKB state, since these constraints 
are linear in momenta, (ii) A Wick rotation has been made to Euclidean 
signature (iii) We have applied a variant of the Baker-Campbell-Hausdorf 
formula in isolation of the energy 'eigenvalue' on the state. 



1 n „ 
n ' i 

(120) 

Here H'^^ is the eigenvalue upon exponentiation of the operator equation 
(112), and is made up of multiple commutators of H with Hd- The precise 
expression is not important so much as the fact that it is highly singular 

oo 

i?^, = J]n'^G"(j(3)(0))X, (121) 

n=0 

due to the factors of ^(^•'(O), which when exponentiated yield zero (for the 
appropriate spacetime signature such that the argument of the exponential 
is negative). Or, if nonzero, oscillates with infinite frequency (for an appro- 
priate signature) such that the contributions from all WKB states cancel. 
Whatever the case, the quantum transition amplitude can be defined so as 
to suppress the WKB states so that they are either highly pathological or 
nonmeasur able . 

Let us assume the appropriate conventions for measuring states, between 
finite time intervals, are in place. Unless the lapse density is zero (either for 
for the time component of the metric to vanish goo = 0, or for the spatial 
3-metric to be degenerate det{hij) = 0), the semiclassical states cannot 
contribute unless r = tf — ti = 0. The problem of time has resurfaced. 
It also implies that the laws of physics for gravity are valid or accessible 
only in the quantum realm, and that we must regard with skepticism the 
semiclassical limit if any sense is to be made of path integral quantization. 
On the other hand, we do live in a universe in which we make contact with 
reality by testing the classical limit of a quantum theory. So what is the 
way out? 

A way out is possibly via the semiclassical-quantum correspondence! In 
other words, when l^'i^jtb) = |^'can)) the path integral is in conformity with 
the canonical approaches to an experimentally verifiable quantum theory of 
gravity. So the question now becomes, for which sets of quantum gravita- 
tional states are the quantum and semiclassical limits identically the same? 
It is the set of generalized Kodama states "^cKod- The canonical quantum 
states are actually a subset of the WKB states [ ], and they are preferen- 
tially selected by the path integral for a Hilbert space upon which to make 
measurements. So one does not necessarily need to access the Planck scale 
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in order to experimentally test quantum gravity. One can do so, vicariously, 
in the semiclassical limit via the SQC. 

The dichotomy between these two extremes exists in the full theory of 
gravity, for which quantum singularities make the split evident. However, 
there are intermediate cases as well, fully consistent with the theory of quan- 
tum gravity. This are precisely the minisuperspace models. 

Minisuperspacc models, albeit an approximation for the full theory of 
quantum gravity (or perhaps vivce versa), are perfectly suitable models in 
their own right modulo issues of the commutativity of superspace reduction 
with quantization. For usch models the S^^^ (0) singularities in the quantum 
terms are no longer present and one may still have for its corresponding 
WKB states relations of the form 

H\^^ini) = Ktl^f^.m) = {HcihGfH + n^G^H,tf^2 + ...) |*™m> / 0. (122) 

These states are deviod of infinities, and performing the analogous insertion 
gives something of the form 

^(^r(Sv),</>(s,,)|^L"L><^Lt.lA"(2*j,</>(s,j) 

n 

n 

= E *^:L[A"(St,), </'(St,)]*2„JA"(s,j, 0(s*j] 

n 

= / DADcPexp{i [ LA.ft[A,</.])(123) 

J JM 

In the case where there is a large set of quantum states of minisuperspace it 
may not be as obvious the relation of a particular one to the path integral 
but at least the left-hand side is finite, which implies that that the path 
integral most likely converges. 



6 Discussion 

The fundamental principle underlying the semiclassical-quantum conjec- 
ture for quantum gravity is that the Hamiltonian constraint which is either 
quadratic or cubic in momenta can be transformed into a form linear in the 
momenta. In this form the semiclassical state can be exactly determined 
by solution of a first-order functional differential equation for the full theory 
without resorting to minisuperspace reduction. If the semiclassical-quantum 
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correspondence holds then the resulting state would also be a quantum state 
and the reduced phase space quantization would be equivalent to the full 
Dirac quantization [13]. Let us briefly compare these two processes as they 
apply to quantum gravity. The reduced phase space quantization, which fac- 
torizes or simplifies the constraint into linear form at the classical level prior 
to promoting the classical operators to quantum form, has the following fea- 
tures: (i) it allows determination of the semiclassical state by solution of a 
first-order functional differential equation, which is easier to solve and can 
generally be solved exactly; (ii) it does not take into consideration quantum 
ambiguities and for this reason is not the exact quantum state. However, 
in the event that the SQC holds, the issue of ordering ambiguities would 
by definition be resolved since any ambiguities would be a-priori isolated in 
the form of quantum counterterms which are dealt with by separately. In a 
certain sense, it is as though the quantum ambiguities are introduced into 
the theory externally by hand through the quantization process. 

The full Dirac quantization procedure, on the other hand, which pro- 
motes the classical operators into their quantum versions while the full con- 
straint is still in unfactorized form, has the following features: (i) it is now 
the quantum state which satisfies a condition, generally a non-linear func- 
tional differential equation. This equation would allow exact determination 
of the quantum state if it could be solved, however due to the infinite number 
of degrees of freedom it is virtually intractable to do so. It is a conjecture 
[8] that the path integral by definition solves the constraint formally for cer- 
tain boundary conditions on the main dynamics variables, however there are 
many properties of of relevance in the path integral which are ill understood, 
for example the transformation properties of the path integral measure, (ii) 
Upon reduction to minisuperspace (which is not the same thing as reduced 
space quantization, and in itself raises separate issues of the communtativ- 
ity of reduction process with respect to quantization), even if it is possible 
to solve the reduced differential equation (which is still non-linear), the or- 
dering ambiguities would still be present, (iii) the ordering ambiguities, 
whether they be in the superspace or in the minisuperspace version of the 
theory, would be reflected in the quantum state and the state would not be 
unique [9] thus even though the range of ambiguity may be parametrizable 
it would still not be clear which quantum theory would be the correct one 
to choose. 

So we see that in general there is a trade-off between both approaches. 

But in particular, the former approach would be preferable provided the 
semiclassical-quantum correspondence holds. It may be possible to relate 
the uncertainties due to the Dirac quantization process, which imposes a 
functional Schrodinger operator equation with ordering ambiguities, to the 
ill-definedness of the path integral measure. It is the natural tendency to 
do so since the integrand of the path integral is classical and involves c- 
numbers. Since there are no operators one may think that any uncertainties 
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can automatically be attributed to the path measure itself which may in 
general contain infinitcfactors.par In Part II and III we will illustrate 
the equivalence of the canonical to path integral approaches to quantization 
when gravity is coupled to matter sources. We will also refine the character- 
ization of the so-called suprious orthonormal states. In [4] we computed the 
antisymmetric contribution to the quantum terms arising from the Hamilto- 
nian constraint. In these next few works we will complete the calculation of 
the remaining terms, and illustrate the algorithm for systematically finding 
the generalized Kodama states correpsonding to the full set of terms. 
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